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Sterile neutrinos in the electronvolt mass range are hinted at by a number of terrestrial neutrino 
experiments. However, such neutrinos are highly incompatible with data from the Cosmic Microwave 
Background and large scale structure. This paper discusses how charging sterile neutrinos under a 
new pseudoscalar interaction can reconcile eV sterile neutrinos with terrestrial neutrino data. We 
show that this model can reconcile eV sterile neutrinos in cosmology, providing a fit to all available 
data which is way better than the standard ACDM model with one additional fully thermalized 
sterile neutrino. In particular it also prefers a value of the Hubble parameter much closer to the 
locally measured value. 
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INTRODUCTION 


Data from the Planck satellite [T and other cosmo¬ 
logical probes have provided astoundingly precise infor¬ 
mation on cosmology. While some cosmological param¬ 
eters, like the spatial curvature, are probed directly by 
the CMB measurements, others can only be inferred in¬ 
directly from a global fit to all cosmological parameters 
simultaneously. The prime example in the last category 
is the Hubble parameter which formally has been deter¬ 
mined very accurately from Planck data [I]. However, 
the inferred value can only be trusted within the stan¬ 
dard ACDM cosmological model. Furthermore there is a 
pronounced discrepancy between the value inferred from 
CMB data and the value measured directly using stan¬ 
dard candles in the nearby universe j5j. This discrepancy 
is mitigated if the model includes additional relativistic 
species (i.e. a non-standard value of N e g). 


Concerning neutrinos, Big Bang Nucleosythesis 
(BBN), CMB and Large Scale Structure (LSS) measure¬ 
ments disfavor a fourth (sterile) neutrino with a mass 
in the eV range, as indicated by oscillation experiments 
m- There are two reasons: On the one hand one addi¬ 
tional degree of freedom is excluded by BBN and CMB 
observations; on the other hand the eV mass scale is ruled 
out with high significance by LSS bounds. Recently, var¬ 
ious models of neutrino self-interactions GHEE! have been 
proposed to reconcile sterile neutrinos in cosmology, by 
preventing additional neutrinos from being fully thermal¬ 
ized in the early Universe. Among the secret interaction 
models, the pseudoscalar model m, where interactions 
are confined to the sterile sector and mediated by a light 
pseudoscalar, can naturally accommodate one eV-sterile 
neutrino in cosmology, solving the tension with BBN and 
CMB measurements. 


MODEL FRAMEWORK AND PRODUCTION OF 
STERILE NEUTRINOS 

Specifically, the model considered couples the mainly 
sterile 4th neutrino mass state, 1 Z 4 , to a new light pseu¬ 
doscalar, <f >, with mass -C leV via 

C ~ g s <f>v/rt 5^4- (1) 

This new interaction provides a background matter po¬ 
tential for neutrinos which is quadratic in g s and, pro¬ 
vided that the coupling constant g s is larger than 10 -6 , 
the sterile neutrino production in the early Universe is 
significantly suppressed until after the decoupling of ac¬ 
tive neutrinos. This means that the total energy density 
in active neutrinos, quantified in terms of 


with po = (4/11 ) 4 // 3 p 7 , can be significantly lower than 
the N e g ~ 4 predicted in the absence of non-standard 
interactions. 

However, even in the presence of a significant non¬ 
standard matter potential, sterile and active neutrinos 
can still equilibrate at later times, and if the sterile neu¬ 
trino has a mass of order 1 eV this can lead to conflict 
with the cosmological bound on neutrino mass [9]. 

Here we will refrain from a detailed calculation of how 
efficient the equilibration is and simply use the most con¬ 
servative assumption possible. Naively we expect that 
full equilibration leads to N e g/A in the sterile sector and 
37Veff/4 in the active sector. However, the pseudoscalar 
interaction leads to a strong coupling between sterile 
neutrinos and pseudoscalars at late times. If this hap¬ 
pens prior to active-sterile equilibration the division will 
be 117Veff/32 in the strongly coupled v s — <j> fluid and 
21iVeff/32 in the remaining non-interacting and massless 
neutrinos. As we will see in the next section structure 
formation data prefer a smaller fraction of the strongly 
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interacting component, so taking the fraction to be 11/32 
is the most conservative choice. 


LATE TIME PHENOMENOLOGY AND FIT TO 
STRUCTURE FORMATION DATA 

The late-time phenomenology of the pseudoscalar 
model leads to a collisional recoupling of sterile neutrinos 
via the process v s v s ut <l><t> before their non-relativistic 
transition (assuming g s > 1CU 6 and m v ^ s > 1 eV) and 
before recombination. In the collisional regime, neutri¬ 
nos and pseudoscalars are not free-streaming, but rather 
behave as a single fluid with no anisotropic stress. The 
impact on the CMB temperature power spectrum con¬ 
sists of an enhancement of the monopole term j|13j , which 
would spoil the Planck measurements if the collisional 
regime was extended to all neutrino species, including 
active neutrinos. 

At temperatures higher than the sterile neutrino mass 
the combined fluid is fully relativistic with an effective 
equation of state parameter of w = 1/3. However, when 
sterile neutrinos go non relativistic, they annihilate into 
pseudoscalars vv —> <fxj>, while the inverse process </></> -A 
vv becomes kinematically prohibited. This increases the 
energy density of the combined fluid relative to that of a 
fully relativistic fluid. During this process the pressure of 
the combined fluid also drops relative to energy density 
because of the importance of the rest mass of the sterile 
neutrino mm- However, for much lower temperatures 
the fluid consists only of pseudoscalars and is again fully 
relativistic with w = 1/3. The redshift where deviations 
in p and w set in depends directly on the sterile neutrino 
mass. We have shown the redshift evolution of p and w 
in the upper and lower panel of Figure^ respectively. 

At the end of the annihilation processes the cosmic 
abundance of sterile neutrinos is highly suppressed and 
there is little suppression of matter fluctuations on small 
scales. This is in stark contrast to the ACDM model 
where the presence of a free streaming massive sterile 
neutrino causes a dramatic suppression of the matter 
power spectrum. 

The difference between the two scenarios can easily be 
seen in Figure [2] where we show the matter power spec¬ 
trum obtained for three different sterile neutrino masses: 
in the pseudoscalar scenario the increase in the sterile 
neutrino mass does not lead to the suppression of power 
at small scales that one finds in the ACDM model. 


ACDM with massive sterile neutrinos 

If the only non-standard physics is the addition of a 
sterile neutrino with the mass and mixing needed to ex¬ 
plain the short baseline data, the expectation is that 
the additional species are almost fully thermalised in the 




FIG. 1: Upper panel: Relative increase in the pseudo scalar- 
sterile neutrino energy density compared to the energy density 
of one active neutrino, due to sterile neutrino annihilations. 
Lower panel: Temporary suppression of the pseudo scalar- 
sterile neutrino equation of state parameter. 

early universe (for a recent treatment see e.g. [23]). The 
prediction therefore is that N e g ~ 4 and m s ~ 1 or 
~ 3 eV [21 E] - Thus, when testing the ACDM model 
with a varying sterile neutrino mass m„ i5 , we will keep 
N e g = 4.046. Hereafter we will refer to this model as 
the A + l^s + m UjS model, while the pure ACDM model 
with 3.046 massless active neutrinos is denoted by A and 
the model with 1 eV fully thermalized sterile neutrino is 
denoted by A + 1 eV v s . 


Pseudoscalar model 

In the pseudoscalar model N e g is a free parameter de¬ 
pending on g s . However, unless the pseudoscalar is ther¬ 
malised by other means in the early Universe the expec¬ 
tation is that 3 < iV e ff < 4 with N e g ~ 3 corresponding 
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FIG. 2: Linear matter power spectra for A + lu B +m„ |B (lower 
panel) and the pseudoscalar model (upper panel) for various 
values of the sterile neutrino mass and for N e g = 4.046 (in the 
pseudoscalar model 11/32 x 4.046 are strongly interacting). 
For comparision, the matter power spectrum obtained with 
Planck 2015 best-fit is also shown. 

t° 9s > 10 -5 and N e g ~ 4 corresponding to g s < 10 -6 . 

Assuming that actives and steriles completely equili¬ 
brate via oscillations prior to the timescales of interest 
to CMB and large scale structure, the expection is that 
11/32 of the total energy density in all neutrinos and 
the pseudoscalar is in the u s — <fi fluid. However, this 
ratio could be smaller if equilibration is incomplete or if 
steriles and pseudoscalars equilibrate after active-sterile 
equilibration. In most runs we fix the ratio, f- ln t, to 11/32, 
but we have also tested the case where /i nt is allowed to 
vary. 

We stress that each value of g s corresponds to one 
value of N e g and one value of but while N e g can 
be found as a function of g s m, the calculation yield¬ 
ing /int is unfeasible to perform using current technol¬ 
ogy, and we simply leave /; n t as a free parameter when 


considering partial equilibration. When also letting N e g 
vary freely, the result is a two-dimensional parameter 
space, where the pseudoscalar models trace out a one- 
dimensional path. 

Comparison with data 

In order to test the ability of the pseudoscalar model 
to fit cosmological data and to constrain its parameters, 
we perform a Markov Chain Monte Carlo analysis using 
CosmoMC m with a modified version of the Boltzman 
solver CAMB D21- Our basic data set is based on Planck 
2015 high multipole temperature data and low multipole 
polarization data (PlanckTT+lowP), implemented ac¬ 
cording to the prescription of Ref. |T5]. Additional data 
sets are the Hubble Space Telescope (HST) prior on the 
Hubble constant from direct measurements [5] and Bary- 
onic Acoustic Oscillations (BAO), including 6dFGS [19] . 
SDSS-MGS HO], BOSS-LOWZ BAO J2T] and CMASS- 

drii [sg. 

X 2 results — 

In order to get a feeling for how well the models fit the 
data compared with each other we have obtained best-fit 
X' 2 values for a variety of different cases (see Table]]]). Be¬ 
fore comparing values for the different models, we should 
state clearly the question we want to ask. We are in 
fact less interested in comparing the pseudoscalar model 
with the standard ACDM model. Instead, we would like 
to compare these two models in a possible future where 
terrestrial experiments have confirmed the existence of 
an eV-scale sterile neutrino. So the two models that 
we should compare are the pseudoscalar model and the 
ACDM model with an additional fully thermalised neu¬ 
trino of mass around 1 eV or higher. 

If only CMB data is used the pseudoscalar model with 
fixed /i n t has a slightly higher x 2 than ACDM and the 
model with varying / int has a slightly lower x 2 - Even 
ACDM with a 1 eV sterile neutrino does not have a sig¬ 
nificantly higher y 2 - However, especially when BAO data 
is included, the sensitivity to the neutrino rest mass in¬ 
creases drastically and x 2 for ACDM with a 1 eV sterile 
neutrino becomes much worse - formally excluding it at 
approximately 5 a relative to pure ACDM. However, the 
pseudoscalar model again has a x 2 comparable to ACDM 
and thus provides a significantly better fit than ACDM 
with a 1 eV sterile neutrino, as shown in Figure [3j In 
conclusion, in the absence of evidence for eV sterile neu¬ 
trinos from short baseline data, cosmological data does 
not merit the inclusion of sterile neutrinos or sterile neu¬ 
trino interactions. However, if the existence of eV sterile 
neutrinos is confirmed by future short baseline data, new 
physics is needed in order to reconcile them with cosmol¬ 
ogy. Our analysis shows that the pseudoscalar interaction 
drastically reduces y 2 and allows for a fit as good as pure 
ACDM. 
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Data 

2 

Xtot 

2 

XCMB 

Xhst 

2 

Xbao 

Model 

CMB 

11265.8 

11265.8 

— 

— 

P 


11259.1 

11259.1 

— 

— 

P + /int 


11274.8 

11274.8 

— 

— 

A + 1 eV v s 


11260.3 

11260.3 

— 

— 

A 

CMB+HST 

11265.8 

11265.8 

0.0 

— 

P 


11260.2 

11260.0 

0.2 

— 

P + /int 


11279.2 

11275.5 

3.7 

— 

A + 1 eV v s 


11266.1 

11262.6 

3.5 

— 

A 

CMB+BAO 

11270.5 

11266.3 

— 

4.2 

p 


11263.8 

11259.4 

— 

4.4 

P + /int 


11288.1 

11279.6 

— 

8.5 

A + 1 eV v a 


11264.7 

11260.2 

— 

4.5 

A 

All 

11272.9 

11267.3 

0.7 

4.9 

P 


11266.0 

11259.6 

1.8 

4.6 

P + /int 


11288.7 

11280.2 

1.7 

6.8 

A + 1 eV v s 


11270.8 

11260.6 

5.8 

4.4 

A 


TABLE I: Best-fit \ 2 f° r various models and various data set combinations. The labels for the various models follow the prescrip¬ 
tion described in the text. Data-set combinations are labeled according to the text. The “All” case refers to CMB+HST+BAO. 
The best-fit x 2 values were obtained with the BOBYQA routine implemented in CosmoMC. The obtained values of \ 2 are typically 
within Ay 2 ~ 1 of the true global best-fit value. 


Data 

n s N eff 

m„, s [eV] I7o[km/s/Mpc] 

/int 

Model 

CMB 

0.978t°;°i* 3.66i°; 2 3 8 6 

Q QCJ + 1 - 1 
O.LM_ 0 .76 

74.0111 

11/32 

P 


0.97791" 3.42!°;“ 

c 7+1.7 
- 2.0 

71.2ll; 5 6 

0 i 1 1+0.050 
u * 111 -0.062 

P + /int 


0.9994 + 0.0052 4.046 

< 0.364 

74 CS + 1 - 6 
'^• O - 0.86 

— 

A + lv s + m„, a 

CMB+HST 

0.977±g;go| g 3.64 + 0.19 

3.00lJ; 2 7 

73.8 + 1.4 

11/32 

p 


0.9848 + 0.0074 3.59 + 0.19 

5.3li;| 

72.8 + 1.4 

0 124+ 0 054 
0.060 

P + /int 


0.9990 ± 0.0045 4.046 

< 0.284 

74.411? 

— 

A + lv s + m v , s 

CMB+BAO 

0-968tg;g“ 7 3.34t°ji 

o fil+0.97 
°- Oi -0.43 

70.75l?;f 

11/32 

P 


0.97061“ < 3-55 

6.0 ±1.9 

69.551 ?? 6 

0 1 04 +0 048 
u - lu ^-0.060 

P + /int 


0.9945 ± 0.0038 4.046 

0.26!°;“ 

72.38 + 0.61 

— 

A + lv s + m„, a 

All 

0.97371 °'0057 3.49 + 0.18 

q 'T'7+0.64 

1 1 -0.51 

71.8 + 1.0 

11/32 

p 


0.9762 ± 0.0058 3.4511g| 

5.711;? 

70.9 + 1.0 

0 12^ +0 051 
u.izo _ 0 057 

P + /int 


0.9946 ± 0.0038 4.046 

0.241°;“ 

72.54 + 0.58 

— 

A + lv s + m„ |S 


TABLE II: Marginalized constraints are given at la, while upper bounds are given at 2 a for the pseudoscalar model (with and 
without varying /j n t) and for the ACDM model with one additional fully thermalized massive neutrino. Data-set combinations 
are the same as in Table |T] 
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FIG. 3: Upper panel: Temperature anisotropies power spec¬ 
tra for various models: Pseudoscalar (red/dashed line), A + 
0.06 eV (blue/solid line), A+lr' a +m„ lS (purple/dot-dashed 
line), A + 1 eV v s (black/solid line). Lower panel: Rel¬ 
ative errors for the models listed above compared to the 
A + 0.06 eV r'a model. In all cases the spectra plotted are 
with the best-fit parameters obtained when fitting to CMB 
data. 


Parameter constraints — 

In Table [TT| the ler marginalized constraints or the 20- 
limits on n s , N e s 7 TO !y . s , Hq are reported for various data 
set combinations and for different models. 

From the analysis it is clear that although N e ff = 4 
is marginally allowed in the pseudoscalar model (within 
2 a when fitting only CMB data), a value below 4 is pre¬ 
ferred. This indicates that relatively high values of g s are 
favored and sterile neutrinos are partially thermalized in 
the early universe. E.g., when BAO are included in the 
analysis and /i„ t is allowed to vary, N e g < 3.55 at 95% 
c.l. This upper bound on N e g translates into a lower 
bound on the coupling g s > 3 x 10 -6 , which is consistent 
with the expectations. 


FIG. 4: One dimensional posterior for „ in the A + lu s + 
s model and in the pseudoscalar scenario for the combi¬ 
nation of data PlanckTT+lowP. The vertical lines show the 
best fit for the sterile neutrino masses obtained by the global 
oscillation data analysis in the 3+1 scenario [2] 

The preference for /; nt > 0 at more than 2er underlines 
that additional neutrino species cannot be free-streaming 
at recoupling, thus, if the short baseline experiments will 
confirm the existence of a sterile neutrino, non standard 
interactions will be required in order to accommodate a 
fourth neutrino in cosmology. 

Note that the usual correlation between the scalar 
spectral index n s and N e g is driven by diffusion damping 
at large l which are absent for a tightly coupled fluid. 
This effect is clearly seen on the lower panel of Figure [3] 
That means that the value of n s is virtually unchanged 
with respect to standard ACDM while the A + 1 eV v s 
model prefers a much higher value as shown in Table |TT]. 

Considering the sterile neutrino mass, when only CMB 
data are considered, m„ iS < 0.364 eV at 2cr, this upper 
limit further tightens when BAO are included; the eV 
range is excluded at high significance in a pure ACDM 
model with free streaming neutrinos. However, the late 
time phenomenology of the pseudoscalar model makes 
sterile neutrinos fully consistent with the eV-mass range: 
the overlap with the results of the global fit in the 3+1 
scenario j2] is not only in the region around Am 2 ~ 1 eV 2 
but also with Am 2 ~ 6eV 2 (see Figure [ 4 ]). 

In Figure [5] we show that models with a strongly self¬ 
interacting dark radiation component provide a remark¬ 
ably good fit to CMB data with a preferred value of the 
Hubble parameter fully consistent with local measure¬ 
ments. 

Finally, Figure [6] shows one-dimensional posteriors and 
1 and 2 er marginalized contours for (N e ff, m yjS , H 0 ), ob¬ 
tained within the pseudoscalar scenario and with various 
data set combinations. 
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FIG. 5: One dimensional posterior for Ho in the ACDM 
model and in the pseudoscalar scenario for the combination 
of PlanckTT+lowP data. The grey region shows the lcr con¬ 
fidence interval from direct measurements. 


DISCUSSION 


We have tested the pseudoscalar model against the 
most precise available cosmological data and found that 
the model is generally compatible with data, providing 
at least as good a fit as the standard ACDM model. Fur¬ 
thermore the fit is vastly better than ACDM with an 
additional sterile neutrino in the eV mass range. 

If the eV sterile neutrino interpretation of short base¬ 
line data turns out to be true cosmology is faced with a 
very serious challenge. Taken at face value such a model 
is excluded by CMB and large scale structure data at 
least at the 5er level. With this in mind it is clear that 
accommodating eV sterile neutrinos requires addition of 
new physics either in cosmology or in the neutrino sector 
(see e.g. [3D] for a discussion). 

The model discussed here provides a simple and ele¬ 
gant way of reconciling eV sterile neutrinos with preci¬ 
sion cosmology. We again stress that this model has a 
late-time phenomenology very different from models with 
purely free-streaming neutrinos and that it could well be 
possible to test details of the model with the greatly en¬ 
hanced precision of future cosmological surveys such as 
Euclid [3T1 . 

Finally, it is interesting that a recent study by Les- 
gourgues et al. [22: find that current cosmological data 
prefers relatively strong self-interactions between dark 
matter and a new dark radiation component. While the 
model presented here cannot provide such dark matter 
interactions at the required strength unless the funda¬ 
mental coupling becomes close to unity, it could be a 
another indication that we are seeing the first signs of 
new, hidden interactions in the dark sector. 


We thank Jan Hamann for comments on the 
manuscript. STH and TT thank the Mainz Institute for 
Theoretical Physics for its hospitality and support dur¬ 
ing the completion of this work. MA acknowledges par¬ 
tial support from the European Union FP7 ITN INVIS¬ 
IBLES (Marie Curie Actions, PITN- GA-2011- 289442) 
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